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Abstract 

The Thr92Ala (rs225014) polymorphism in the type 2 deiodinase (D/02) gene has been associated with insulin resistance (IR) 
and decreased enzyme activity in human tissues but kinetic studies failed to detect changes in the mutant enzyme, 
suggesting that this variant might be a marker of abnormal DI02 expression. Thus, we aimed to investigate whether other 
DI02 polymorphisms, individually or in combination with the Thr92Ala, may contribute to IR. The entire coding-region of 
DI02 gene was sequenced in 12 patients with type 2 diabetes mellitus (T2DM). Potentially informative variants were 
evaluated in 1077 T2DM patients and 516 nondiabetic subjects. IR was evaluated using the homeostasis model assessment 
(HOMA-IR) index. DI02 gene sequencing revealed no new mutation but 5 previously described single nucleotide 
polymorphisms (SNPs). We observed that all T2DM patients displaying high HOMA-IR index (n = 6) were homozygous for the 
rs225017 (T/A) polymorphism. Further analysis showed that the median fasting plasma insulin and HOMA-IR of T2DM 
patients carrying the TAT genotype were higher than in patients carrying the A allele (P = 0.013 and P = 0.002, respectively). 
These associations were magnified in the presence of the Ala92Ala genotype of the Thr92Ala polymorphism. Moreover, the 
rs22501 7 and the Thr92Ala polymorphisms were in partial linkage disequilibrium (|D'| = 0.81 1; ^ = 0.365). In conclusion, the 
rs225017 polymorphism is associated with greater IR in T2DM and it seems to interact with the Thr92Ala polymorphism in 
the modulation of IR. 
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Introduction 

Thyroid hormones are critical to the development, growth and 
metabolism of virtually all tissues [1,2]. The iodothyronine 
deiodinases types 1 (Dl), 2 (D2) and 3 (D3) are selenoenzymes 
of the oxidoreductase family which catalyze iodine removal from 
the outer (Dl and D2) or inner (Dl and D3) ring of thyroid 
hormones. While Dl and D2 convert T4 to the metabolically 
active hormone T3, D3 inactivates both T4 and T3 [2]. In 
humans, D2 is the most important deiodinase for intracellular T3 
generation in target tissues [1,2] and, together with Dl, 
contributes for 80% of peripheral T3 [3]. The gene that encodes 
D2 (DI02) is expressed in the thyroid, pituitary, brain, heart, 
placenta, skeletal muscle and adipocytes [4,5,6,7,8,9] . D2 activity 
is tightly regulated at transcriptional, post-transcriptional and post- 
translational levels [10,1 1], which thereby supports the hypothesis 
that this enzyme plays an important homeostatic role in 
metabolism. 

Type 2 diabetes mellitus (T2DM) is a heterogeneous group of 
disorders usually characterized by varying degrees of insulin 
insufficiency and insulin resistance (IR), which result in increased 
blood glucose concentrations. Ultimately, IR results either from 



inappropriately increased hepatic gluconeogenesis or decreased 
glucose disposal rate in tissues such as skeletal muscle, adipose 
tissue and liver [12]. 

Glucose homeostasis may be affected by thyroid status [13,14]. 
Thyroid hormone affects insulin action in skeletal muscle and 
adipose tissue in part by upregulating the expression of the glucose 
transporter 4 (GLUT4) and, therefore, increasing glucose uptake 
[15,16]. An increase in both hepatic endogenous glucose disposal 
and insulin mediated-glucose uptake is observed in patients with 
hyperthyroidism when compared with euthyroid subjects 
[17,18,19]. In contrast, studies showed that, in hypothyroidism, 
there is a decrease in both glucose disposal [19] and insulin 
mediated-glucose uptake in muscle, and also an impaired ability of 
insulin to stimulate glucose disposal related to insulinaemia 
[19,20]. In animal models, experimental induction or spontaneous 
forms of thyroid dysfunction are also associated with impaired 
glucose tolerance [20] . Mice with targeted disruption of the dio2 
gene have higher insulin levels during glucose tolerance tests and 
reduced glucose uptake during insulin tolerance tests, consistent 
with the occurrence of IR [21]. According to these findings, it is 
plausible to postulate that a lower intracellular D2-generated T3 in 
skeletal muscle might create a state of relative intracellular 
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hypothyroidism, decreasing the expression of genes involved in 
energy use, such as GLUT4, and thus resulting in increased IR. 

Previous studies have demonstrated that polymorphisms in the 
DI02 gene might interfere in the phenotypic expression of D2 
[7,22]. The DI02 single-nucleotide polymorphism (SNP), in 
which a threonine (Thr) changes to alanine (Ala) at codon 92 
(Thr92Ala, rs2250 14, A/G), was associated with IR in obese white 
women and with a 20% lower glucose disposal rate in white non- 
diabetic women [23]. The Ala92Ala genotype was also associated 
with a more pronounced IR in patients with T2DM and the 
frequency of this variant allele was found to be increased in some 
ethnic groups, such as Pima Indians and Mexican-Americans, who 
have a higher prevalence of IR [23]. Even though some studies 
have failed to demonstrate an association between the Thr92Ala 
polymorphism and glycemic traits or T2DM [24,25,26], the 
association between this polymorphism and T2DM is supported 
by data reported in a recent systematic review and meta-analysis 
[27]. However, the mechanism by which this occurs is still not 
clear. Decreased D2 activity has been observed in human tissue 
biopsy samples [7,27] but studies failed to identify any changes in 
the ex vivo biochemical properties of the mutant enzyme [7,28], 
suggesting that this variant might be only a marker for abnormal 
DI02 expression. Here, we aimed to determine whether other 
SNPs in the DI02 gene, individually or in combination with the 
Thr92Ala polymorphism, contribute to the IR phenotype. 

Materials and Methods 

Ethics statement 

The information obtained during the study did not affect the 
patients' diagnosis or treatment. The protocol was approved by 
the Committee on Research Ethics from Hospital de Clinicas de 
Porto Alegre, and all subjects signed an informed consent term. 
Clinical investigation was conducted according to the principles 
expressed in the Declaration of Helsinki. 

Study population 

In the first phase of this study, we used direct automated 
sequencing of the DI02 gene to search for SNPs in 1 2 T2DM 
patients with different degrees of IR. These 12 patients were 
selected from two subgroups in our T2DM population: 6 patients 
with a low HOMA-IR index (ranging from 0.32 to 1.3) and 6 
patients with a high HOMA-IR index (6.2-23.4). The results 
generated by sequencing the DI02 gene in T2DM patients were 
compared with the sequences of the human DI02 gene available 
in GenBank (http://www.ncbi.nlm.nih.gov/genbank/). 

In the second phase of the study, allele and genotype 
frequencies of a potentially informative DI02 variants found by 
sequencing were compared between 1077 patients with T2DM 
and 516 nondiabetic subjects. The criteria used to select the 
variant for further analysis were: 1) the polymorphism was not 
known to be a neutral variation; 2) it involved change between 
chemically different amino acids or was localized at regulatory 
regions, such as the promoter or 3'- untranslated region (3'UTR); 
3) its frequency was not low (minor allele frequency less than 5% in 
the literature or GenBank); and 4) the frequency of the variant was 
different in the groups under analysis (in this case, low and high 
HOMA-IR groups) [29,30]. Additionally, haplotypes constructed 
from the combination of this polymorphism with the Thr92Ala 
polymorphism were analyzed regarding their effect on IR in 
T2DM patients. 

The 1077 patients with T2DM were selected from a multicenter 
study that started recruiting patients in southern Brazil in 2002 
and whose purpose was to investigate risk factors for T2DM and 



its complications. Initially, it included 4 tertiary-teaching hospitals 
in the Brazilian state of Rio Grande do Sul, namely Grupo 
Hospitalar Conceicao, Hospital Sao Vicente de Paula, Hospital 
Universitario de Rio Grande and Hospital de Clinicas de Porto 
Alegre. The detailed description of that study can be found 
elsewhere [31]. The group of subjects without diabetes was 
composed of 516 healthy volunteers who came to Hospital de 
Clinicas de Porto Alegre to donate blood. 

Clinical and anthropometric profiles and laboratory 
analyses 

A standard questionnaire was used to collect information from 
all patients about age, age at T2DM diagnosis and drug treatment. 
All patients with T2DM were white. The ethnic group was defined 
on the basis of self-classification and subjective classification (skin 
color, nose and lip shapes, hair texture and family history). 

All T2DM patients underwent physical and laboratory evalu- 
ations. They were weighed (barefoot and wearing light outdoor 
clothing) and had their height measured. Body mass index (BMI) 
was calculated as weight (kg) divided by height squared (m 2 ). Blood 
pressure (BP) was measured twice, in the sitting position, with a 5- 
min rest between measurements, using a mercury sphygmoma- 
nometer (Korotkolf phases I and V). The mean of the two 
measurements was used to record systolic and diastolic BP. 
Arterial hypertension was defined as BP& 140/90 mmHg or use of 
antihypertensive drugs regardless of BP at the time of assessment. 

Serum samples for laboratory tests were collected after a 12-h 
fast from all patients with T2DM. Glucose levels were determined 
by a glucose oxidase method; creatinine by the Jaffe reaction; 
glycated haemoglobin (A1C) by an ion-exchange HPLC proce- 
dure (Merck-Hitachi L-9100 glycated haemoglobin analyzer, 
Merck, Darmstadt, Germany; reference range: 2.7-6.0%); and 
total plasma cholesterol, HDL-cholesterol, and triglycerides, by 
enzymatic methods. LDL-cholesterol was calculated using the 
Friedewald equation. Serum insulin was measured by radioassay 
(Elecsys Systems 1010/2010/modular analytics El 70, Roche 
Diagnostics, Indianapolis, IN). Insulin sensitivity was estimated 
by the HOMA-IR index [fasting insulin (milliunits per millilitre) x 
fasting glucose (millimoles per litre)/22.5] [32]. The mean 
HOMA-IR index of control subjects in our laboratory was 
1.84±1.02 [33]. 

Molecular analysis 

DNA was extracted from peripheral blood leukocytes by a 
standardized salting-out procedure. A search for variants of the 
DI02 gene was performed in 12 patients with T2DM by direct 
sequencing of all exons, partial intron sequences (500 bp of exon- 
intron junctions), and 5'UTR and 3'UTR sequences (1000 bp 
each, including the selenocysteine insertion sequence element - 
ESECIS - in the 3'UTR). Screening of this limited number of 
individuals may fail to detect some rarer polymorphisms; however, 
this number seems to be adequate to identify representative 
variants which are sulficiently polymorphic to warrant association 
studies. Direct sequencing in an automated ABI 3100 Avant 
Genetic Analyzer (Life Technologies, Foster City, CA, USA) was 
performed using ABI Prism Big Dye Terminator Cycle Sequence 
Ready reaction kit (Life Technologies) according to the manufac- 
turers' recommendations, and using primers described in Table 
SI. 

Potentially informative variants identified through the sequenc- 
ing of the DI02 gene were selected for subsequent genotyping in 
all patients with T2DM and nondiabetic subjects. These SNPs, 
together with the rs225014 (Thr92Ala) polymorphism, were 
determined using primers and probes contained in Human 
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Custom TaqMan Genotyping Assays 40x (Assays-By-Design 
Service, Life Technologies). Sequences of primers and probes 
used for genotyping are shown in Table SI. Reactions were 
conducted in 96-well plates, at a total reaction volume of 5 U.L and 
using 2 ng of genomic DNA, TaqMan Genotyping Master Mix 1 x 
(Life Technologies), and Custom TaqMan Genotyping Assay lx. 
Plates were then placed in a real-time PGR thermal cycler (7500 
Fast Real-Time PCR System, Life Technologies) and heated for 
10 minutes at 95°C, followed by 45 cycles of 95°C for 15 seconds 
and 62°C for 1 minute. Fluorescence data files of each plate were 
analyzed using automated allele-calling software (SDS 2.1, Life 
Technologies). 

Statistical analyses 

Results are expressed as mean ± SD, % or median (minimum- 
maximum values). Allelic frequencies were determined by gene 
counting, and departures from the Hardy-Weinberg equilibrium 
(HWE) were investigated using chi-squared test. Linkage disequi- 
librium (LD) between the candidate polymorphism and the 
Thr92Ala polymorphism was examined using Lewontin's D' 
|D'| and r measures [34,35]. Haplotypes constructed from the 
combination of these DI02 polymorphisms and their frequencies 
were inferred using the Phase 2.1 program, which uses a Bayesian 
statistical method [35,36]. We also used this program to compare 
the distributions oiDI02 haplotypes between patients with T2DM 
and nondiabetic subjects by permutation analyses of 1,000 random 
replicates [35]. 

Clinical and laboratory characteristics were compared using 
ANOVA, unpaired Student t test or chi-squared test, as 
appropriate. Variables with a skewed distribution were logarith- 
mically transformed before analyses. To check if the candidate 
polymorphism was independently associated with T2DM, a 
multiple logistic regression analysis was performed with T2DM 
as the dependent variable and age, sex, and the possible 
informative polymorphism as independent variables. A linear 
regression analysis was performed to evaluate the independent 
association between the candidate polymorphism and IR param- 
eters after adjusting for covariates (age, sex, BMI and use of 
T2DM medication). Haplotype interaction between the candidate 
polymorphism and the Thr92Ala polymorphism in modulating 
fasting insulin and HOMA-IR index was tested by general linear 
model univariate analyses (GLM), after adjusting for covariates 
(age, sex, BMI and use of T2DM medication). In these interaction 
analyses, each of the polymorphisms was modeled as a dichoto- 
mous variable. 

All analyses were performed using the SPSS 18.0 (SPSS, 
Chicago, IL, USA). The level of statistical significance was set at P 
less than 0.05. Power calculations for the case-control study were 
done using the software PEPI, version 4.0, and showed that this 
study has a power of approximately 80% at a significance level of 
0.05 to detect an OR of 1 .35 or higher (for a 40% frequency of the 
minor allele). 

Results 

DI02 screening and identification of candidate 
polymorphisms for association with insulin resistance 

An interim study comprising samples of 1 2 patients with T2DM 
was used to search for polymorphisms in the DI02 gene. This first 
phase study identified 448 amplicons (data not shown). Although 
no new mutation was detected in the data set, 5 previously 
described SNPs were identified (Figure 1). These SNPs included 
two polymorphisms in the 5' flanking region (rsl99598135, 
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Patients Sax Age Candidate Polymorphisms HON4A-IR index 

rsl9959S135 (-/T] rsl2S85300 (Gly/Asp)- rs225014 (Thr/Ala) rs225015 (A/G) CS225017 (A/T) 

High HOMA-IR index 

41 M 55 T Gly Als A T 23.40 

«2 M 71 Gly Ada AT 2009 

43 M 45 - Gly Ala G T 19.07 

44 F 52 - Asp Ala A T 7.78 

45 M 38 - Gly Thr A T 7.54 

46 F 66 Gly Ala G T 6.20 
Frequency (If) - - 1/6 (16) 5/E [S3 ) 5/6 (83 ) 4/6 167 ) 6/6 (100) 



#7 F 48 - Gly Ala A T 1.33 

48 M 64 Asp Thr A A 1.22 

410 M 59 - Asp Thr G A 0.48 

#12 M 49 - Asp Thr A A 0.32 

Frequency PS) ■ 0/6 (0) 3/6 |50) 2/6 (33) 4/6 (67) 2/6 (33) 



Figure 1. Candidate polymorphisms identified by sequencing 
of the DI02 gene. A. The vertical arrows show the five candidate 
variants in human DI02 gene identified through sequencing. Black 
boxes are coding regions, 'polymorphisms associated with T2DM/IR/ 
fasting insulin. ESECIS, Selenocysteine Insertion Sequence Element; 
3' UTR = 3 '-untranslated region. B. The characteristics of the 12 patients 
with T2DM selected for the screening of the DI02 gene. These patients 
had extreme HOMA-IR indexes and were selected from two subgroups 
according low or high HOMA-IR values. 
doi:1 0.1 371/journal.pone.01 03960.g001 



rsl 2885300), one polymorphism in exon 3 (rs225014 - Thr92/\fa), 
and two in the 3'UTR (rs225015, rs225017). 

The characteristics of the 1 2 patients with T2DM selected for 
the screening of the DI02 gene, as well as the occurrence of the 5 
DI02 polymorphisms identified in these patients are shown in 
Figure 1. We observed that the T/T genotype of the rs225017 
(A/T) polymorphism was more frequent in patients with a high 
HOMA-IR index (6 in 6) than in patients with a low HOMA-IR 
index (2 in 6). Moreover, the relative frequency of this 
polymorphism for white subjects in GenBank was 0.36 whereas 
in our sample was 0.66 (P = 0.129). This polymorphism also fills 
the criteria for selecting an informative polymorphism (see 
Material and Methods Section) and was then selected for further 
analysis in an association study. The Thr92Ala polymorphism also 
fulfilled the above criteria; however, this polymorphism was 
already reported in our sample population as being associated with 
both T2DM risk [27] and increased IR and insulin levels [37]. 
Thus, in the present study, we aimed to evaluate whether 
haplotypes constructed from the combination of Thr92jAa 
polymorphism with the rs225017 polymorphism are associated 
with IR in T2DM patients. 

The rs225017 variant and T2DM development in a case- 
control study 

The baseline characteristics of the 1077 patients with T2DM 
and the 516 nondiabetic subjects included in the study were as 
follows: mean age was 59. 2± 10.1 years vs. 46. 2± 10.3 years (P< 
0.001), and males comprised 47% and 63% of the sample (P< 
0.001). 

The T allele frequency of the rs225017 polymorphism did not 
differ significantly between patients with T2DM and nondiabetic 
subjects (0.51 vs. 0.48, P = 0.196). In addition, genotypes 
frequencies of this polymorphism were similar between T2DM 
patients and nondiabetic subjects (25.9% A/ A, 47.2% A/T, 
26.9% T/T vs. 28.5% A/ A, 46.9% A/T, 24.6% T/T; P = 0.452). 
These findings did not change after logistic regression analysis 
adjusting for age and sex (OR=1.15, 95% CI = 0.86-1.55; 
P = 0.354). Genotype frequencies of the rs225017 polymorphism 
were in HWE in patients with T2DM and nondiabetic subjects 
(P>0.05). 
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Table 1. Clinical and laboratory characteristics of patients with T2DM broken down according to D/02 rs225017 (T/A) 
polymorphism. 



Genotypes 




Total (n = 1077) 


A/A or A/T (n = 798) 


T/T (n = 279) 


P 


r\ge ^rearsj 


59 2 + 1 0 1 


59 6 + 9 8 


58 2 -1- 1 0 3 


0 051 


Duration of diabetes (Years) 


1 1 7 + 8 8 


12i+go 


10 5 + 7 9 


0 009 


Sex (% males) 


47 0 


49 0 


42 0 


0 038 


Dlvll ^Ky/m ) 


28 9 + 4 9 


28 8 + 4 8 


29 3~*~5 3 


0 117 


Total cholesterol (mmol/l) 


5 4+12 


5 4+12 


5 5 + 12 


0 419 


jysionc or irnrn nyj 


1 41 g + 23 2 


1 41 9 + 23 7 


I *f I .o_z I .o 


0 953 


Diastolic BP (mm Hg) 


85.0±12.8 


84.8 ±12.7 


85.6±13.1 


0.823 


HDL cholesterol (mmol/l) 


1.17±0.3 


1.17±0.3 


1.16±0,3 


0.753 


LDL cholesterol (mmol/l) 


1.29±1.1 


1.29±1.1 


1 .31 ±1.1 


0.533 


Serum creatinine (fimol/l) a 


102.8 (35.4-937) 


81.4 (26.5-937) 


79.6 (35.4-540) 


0.127 


Triglycerides (mmol/l) a 


2.00 (0.3-13.2) 


2.00 (0.3-13.2) 


1 .99 (0.5-8.7) 


0.821 


Fasting plasma glucose (mmol/l) 


9.4 (2.9-21.9) 


9.3 (2.9-21.9) 


9.8 (3.3-20.5) 


0.144 


SU/Met/SU + Met (%) 


32.5/38.5/16.5 


31.1/36.0/14.3 


34.8/45.9/22.9 


0.249/0.027/0.015 


A1C (%) 


7.2±2.0 


7.1 ±2.0 


7.3 ±1.8 


0.343 


Fasting insulin (Ul/ml) a ' b 


13.9 (1-46.1) 


10.6 (1-44.9) 


15.7 (3.2-46.1) 


0.005 


HOMA-IR index a ' b 


3.8 (0.3-20.1) 


3.5 (0.3-19.5) 


5.2 (0.6-20.1) 


0.005 



Data are mean ± SD, median (minimum - maximum values) or %. A1C, glycated haemoglobin; BMI, body mass index; HOMA-IR, homeostasis model assessment - insulin 
resistance; Met, metformin; SU, sulfonylureas; WHR, waist-to-hip ratio; P values were estimated by % 2 or ANOVA, as appropriate. 
a Variables which were logarithmically transformed before analyses. 

b For comparisons of fasting insulin levels and HOMA-IR index among rs22501 7 genotypes, we analyzed only 227 individuals (162 individuals harboring the A/A or A/T 
genotypes and 65 individuals harboring the T/T genotype). 
doi:1 0.1 371 /journal.pone.01 03960.t001 



The rs225017 DI02 polymorphism is associated with 
insulin resistance in patients with T2DM 

Table 1 summarizes the clinical and laboratory characteristics 
of patients with T2DM grouped according to the different 
genotypes of the rs225017 polymorphism. Assuming a recessive 
model of inheritance, patients with A/ A or A/T genotypes were 
grouped and compared with patients carrying the T/T genotype. 
T2DM duration (P = 0.009), sex proportion (P = 0.038) and use of 
metformin alone (P = 0.027) or metformin + sulfonylurea 
(P = 0.015) for treatment of T2DM were differentially distributed 
between patients with the T/T genotype and A allele carriers. It is 
worth mentioning that none of the variables showed in Table 1 
attained statistical significance when assuming dominant (A/ A vs. 
A/T-T/T) or additive (A/A vs. T/T) models of inheritance (data 
not shown). 

A subgroup of 227 patients with T2DM who were not receiving 
insulin was further analyzed for IR measurements. This subgroup 
was representative of the whole sample: mean age was 59.3±9.3 
years (P = 0.803 for the comparison with the whole sample), mean 
T2DM duration was 11.2±8.2 years (P = 0.580), mean A1C was 
7.3±1.7% (P=0.12), and mean BMI was 28.9±4.7 kg/m 2 
(P = 0.89). Males comprised 45.3% (n= 103) of the sample 
(P = 0.45). 

In this subgroup of patients, median fasting plasma insulin levels 
were higher in patients with the T/T genotype than in patients 
carrying the A allele [15.71 (3.2-46.1) vs. 10.57 (1-44.9), 
P = 0.005], whereas fasting plasma glucose levels did not differ 
significantly between groups (P = 0. 144). Moreover, T/T genotype 
patients showed a higher HOMA-IR index than patients carrying 
the A allele [5.20 (0.6-20.1) vs. 3.50 (0.3-19.5), P = 0.005]. 
Because insulin sensitivity is known to be affected by multiple 



independent factors, multiple linear regression analyses were 
performed with HOMA-IR index (login HOMA-IR) or insulin 
level (log 10 fasting insulin) as dependent variables. The T/T 
genotype of the rs225017 polymorphism remained significandy 
associated with HOMA-IR [standardized coefficient B for T/T 
genotype =-0.366, 95% CI (-0.654—0.078); P = 0.013] and 
insulin levels [standardized coefficient B for T/T genotype = — 
0.268, 95% CI (-0.441—0.096); P = 0.002], after adjusting for 
age, sex, BMI and use of medication for T2DM (mainly 
metformin or sulfonylureas). 

Association of the DI02 rs225017 and Thr92Ala variants 
with insulin resistance 

We used a Bayesian statistical method to estimate the 
frequencies of different haplotypes produced by the combination 
of the rs225017 and Thr92Ala polymorphisms in patients with 
T2DM. The first letter of the haplotypes refers to the Thr92Ala 
polymorphism, and the second, to the rs225017 polymorphism. 
The wild type haplotype (Thr/A) was observed in 49.0% of the 
sample, the Thr/T haplotype was observed in 18.0%, the Ala/A 
haplotype in 4.0%, and the Ala/T haplotype (double mutated) in 
29.0% of the sample. 

Taking these results into account, we tested whether any 
haplotype constituted by the rs225017 and Thr92Ala polymor- 
phisms might affect fasting plasma insulin levels or HOMA-IR 
index differently as compared with the effect of each polymor- 
phism analyzed separately. To test this hypothesis, we performed 
GLM analyses using fasting insulin levels or the HOMA-IR index 
as the dependent variable, and age, sex, BMI, use of medication 
for T2DM and the polymorphic combinations of the two DI02 
variants as independent variables. These interaction analyses are 



PLOS ONE | www.plosone.org 



4 



August 2014 | Volume 9 | Issue 8 | e103960 



New Polymorphism in DI02 Is Associated with Insulin Resistance 




AAor AT 



rs225014 (Thr92Ala) 

I I Thr-AlaorThr-Thr 
Q Ala-Ala 



rs225017 (A/T) 



00 CTi 



rN ro i- 



Figure 2. Possible interaction between rs225017 and Thr92Ala 
variants in modulation of insulin resistance. HOMA-IR indexes in 
T2DM patients according to different genotypes of DI02 rs225017 (A/T) 
and rs225014 (Thr92Ala) polymorphisms. Results are expressed as 
median (percentiles 25% and 75%). Circles represent outlier values. 
Interaction P value =0.010, adjusted for age, sex, BMI, and use of 
medication for T2DM. 
doklO.1 371/journal.pone.01 03960.g002 
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shown in Table 2. Patients carrying the Ala/ Ala - T/T 
haplotype showed higher fasting insulin values than patients with 
other genotype combinations, adjusting for covariables 
(F= 11.072, P = 0.001). Likewise, patients carrying the Ala/ Ala - 
T/T haplotype showed a higher HOMA-IR index than 
patients with other genotype combinations (F = 4.740; P = 0.010) 
(Figure 2). 

To determine whether the rs225017 (A/T) and rs225014 
(Thr92Ala) polymorphisms were in LD in T2DM patients, we 
used Lewontin's |D' | and r z measures. The analysis showed that 
the rs225017 polymorphism was in partial LD with the rs225014 
polymorphism (|D'| =0.811; r 2 = 0.365). These results were 
confirmed in the HAPMAP Project (http://hapmap.ncbi.nlm. 
nih.gov/). 



x 
< 



II 



r- <J\ m 



("N CN CO 



Oi Oi » ^ 



E % 



9 o £ 



Q. ">_ 
Q. -7 

« a 



3 QT JJ 



■n ° 



E 

=3 r 

.1 - - - 

x o u £ 

to „ c h 

F « = £■ 

aj — 

i *j Ji; 

= *- a) 

i E £ 

.l-S » 

c d aj 

*E *5=» ■£ 

S£s 

m c tj 



fin'" 

g ° 3 

= is o 

,y a cc d 

C T ifl 



: m 



O i- ^ QJ 
— .One 



+i s. ° 



5 -o -a q- 

(T3 (TJ C 



S CC CC is 



■a ^ 



D P J 



nj c .£5 3 3 



73 "° m 
CU OJ " 



.2. 2 
o 



' 5 is 



Discussion 

Insulin resistance is a heterogeneous disorder influenced by 
multiple environmental and genetic factors [38,39,40]. Previous 
studies indicate that the Thr92Ala polymorphism in the DI02 
gene is involved in IR pathogenesis [7,23,37]. Even though the 
Ala/ Ala genotype was associated with lower D2 activity in thyroid 
biopsy samples, in vitro studies did not show any change in the 
biochemical properties of the mutated enzyme, suggesting that this 
SNP might be only a polymorphic marker [7]. Therefore, here we 
searched for other informative SNPs in the DI02 gene which 
might be associated with IR or T2DM. We report a new DI02 
polymorphism associated with increased IR in T2DM patients. 
Homozygosis for the allele T of the rs225017 (A/T) polymorphism 
was associated with IR in patients with T2DM, whereas patients 
carrying a combination of the T/T genotype of the rs225017 
polymorphism and the Ala/Ala genotype of the Thr92Ala 
polymorphism showed an increased HOMA-IR index as com- 
pared with patients with other genotype combinations, suggesting 
that they might interact in the modulation of IR. 

The rs205017 (A/T) polymorphism is located at the 3'UTR of 
the DI02 gene. In selenoenzymes, the 3'UTR contains the 
selenocysteine insertion sequence (SECIS), which is very important 
for looping formation in the mechanism of cis-acting mRNA 
structure of selenocysteine insertion [41,42]. The 3'UTR is also 
the main site for ligation of microRNAs (miRNAs) that can 
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function as translational repressors or triggers of transcript 
degradation by partially pairing to this region in target mRNAs 
[43]. Currently, it is universally recognized that miRNAs are 
major regulators of gene expression and key controllers of several 
biologic and pathologic processes, as such those involved in the 
T2DM pathogenesis [44] . Hence, it is reasonable to speculate that 
the rs205017 polymorphism could be involved in a ligation site for 
a given miRNA, decreasing the DI02 expression. Therefore, we 
used the TargetScan software to look for potential miRNA targets 
in the DI02 gene [45]. Six predicted interaction regions with 
miRNA families were found in the DI02 sequence: miR-9, miR- 
29abcd, miR-30abcdef, miR-193, miR-203 and miR-216a. 
Interestingly, miR-9, miR-29 and miR-30 seem to have a role in 
insulin resistance [44,46]. Although the rs225017 SNP is not 
located at any of these miRNA regions, it is close to some of them 
and, therefore, we could not rule out a potential effect of this SNP 
on DI02 gene expression. Alternatively, it could be obstructing the 
mechanism of selenocysteine insertion and, thereby, reducing the 
amount of D2 available for the activation of T3. In this context, 
polymorphic D2, generated from the Ala/T haplotype sequence, 
might generate less T3 in skeletal muscle and adipocytes, which, 
consequendy, may create a local state of intracellular hypothy- 
roidism, decreasing the expression of genes involved in metabo- 
lism, such as GLUT4, and leading to IR. 

Recently, the DI02 rs65 74549 polymorphism, also located at 
the 3'UTR, was associated with elevated fasting insulin levels and 
higher insulin action, but not with T2DM in Pima Indians [47]. 
This polymorphism may also be a candidate for the functional 
variant associated with IR. The rs6574549 polymorphism is in 
strong LD with the Thr92Ala [47] and rs205017 [48] polymor- 
phisms (LD>0.8; HAPMAP Project, http://hapmap.ncbi.nlm. 
nih.gov/). Therefore, further studies are necessary to determine: 1) 
which polymorphism is the functional variant interacting with the 
Thr92Ala polymorphism in the modulation of IR; 2) whether 
there is another functional polymorphism in the DI02 gene yet to 
be identified; and 3) whether these polymorphisms have combined 
effects on D2 function. The results presented here add to the 
understanding of the molecular interactions controlling IR and 
might partially explain the association with the Thr92Ala 
polymorphism. 

Some factors unrelated to the rs205017 and Thr92Ala 
polymorphisms may have interfered with the findings of the 
present study. First, medications for T2DM treatment may have 
played a role because some are known to affect insulin sensitivity. 
However, we minimized such possibility by excluding insulin- 
treated patients from the group of 227 patients analyzed for IR. 
Furthermore, the rs205017 variant and the Ala/Ala-T/T haplo- 
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